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ABSTRACT 

We present a comprehensive treatment of the spectrum of electric dipole emission 
from spinning dust grains, updating the commonly used model of Draine & Lazarian. 
Grain angular velocity distributions are computed using the Fokkcr-Planck equation; 
we revisit the drift and diffusion coefficients for the major torques on the grain, in- 
cluding collisions, grain-plasma interactions, and infrared emission. We use updated 
grain optical properties and size distributions. The theoretical formalism is imple- 
mented in the companion code, SpDust, which is publicly available. The effect of 
some environmental and grain parameters on the emissivity is shown and analysed. 
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1 INTRODUCTION 

Observational cosmology has entered an area of high pre- 
cision, exemplified by the most recent temperature re- 
sults from se nsitive cosmic microwave background (CMB) 
experiments ([Dickinson et al.l 2004 : IReadhead et al.l |2004| : 



IKuo et al.ll2007l : iHinshaw et aLlbOOSl ). However, foreground 
separation and removal r emains a major chal l enge for any 
CMB measurement (e.g. lEriksen et al.l l2008l : iLeach et al.l 

|2008| ). In addition to the standard Galactic foregrounds, 
free-free, synchrotron and thermal dust emission, an un- 
known "anomalous" dust-correlated emission has been ob- 
served over the last decade, in the microwave region of the 
spectrum. The anomalous emissions was firs t interpreted 
, as free -free e mission fro m shock-heated gas bv lLeitch et al.l 
but iDraine &: L azarian (1998a) showed that this 
would require an extremely high plasma temperature and 
a corresponding unrealistic energy injection rate. They 
proposed instead two possible mechanisms to explain the 
anomalous microwave emission. One of them is the magnetic 
dipole emission from thermal fiu ctuations in the magnetisa - 
tion of interstellar dust grains (|Draine fc LazarianI Il999l ). 
The other possible mechanism, on which the present work 
focuses, is electric dipole radi ation from the sma ll est car - 
bonaceous grains, described in iDraine fc LazarianI (|l998bD . 
hereafter DL98b. The physical principle is quite straight- 
forward: dust grains are presumably asymmetric, and thus 
will have a nonzero electric dipole moment. These grains 
will spin due to interaction with the ambient interstellar 
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medium (ISM) and radiation field, and thus radiate electro- 
magnetic waves due to the rotation of their electric dipole 
moment. To get the electric dipole radiation spectrum, one 
thus needs three ingredients: the quantity of small grains, 
then their dipole moment, and finally their rotation rates. 

Although the observational interest in electric dipole 
radiation from spinning dust grains has only grown in the 
last decade, there is a long standing history of theoretical 
work on the subject. lEricksorJ (|l957h was the first to con- 
sider the possibility that rotatin g dust grains could be the 
source of non-thermal radio-noise. iHovle fc Wickr amasinghd 
(1970) showed that this process was dominated by grains 
with radius a < 10~^ cm and could l ead to radio emission 
around 10 GHz. iFerrara fc DettmaJ (| 19941 ) estimated the 
spinning dust emissivity for thermally rotating grains. The 
first to provide a detail ed treatm en t of ro tational excitation 
of small grains were iRouan et al] (|l992l ') . They considered 
the effect of collision s with gas atoms and absorption and 
emission of radiation. I Anderson fc Watson! (|l993l ) evaluated 
the effect of collisions with ions and "plasma drag" (torques 
due to the electric field of passing ions). 

DL98b provided the first comprehensive study of the 
rotational dynamics of small grains, including all the pre- 
vious effects. They evaluated, as a function of grain radius 
and environmental conditions, rotational damping and ex- 
citation rates through collisions, "plasma drag", infrared 
emission, emission of electric dipole radiation, photoelec- 
tric emission and formation of H2 molecules. The spectra 
they provided are now widely used in interpreting ISM mi 
crowave emission (e.g . iFinkbeiner ,2004: Watson et al. 2005 
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(e.g. iBandav et al. 
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l2007l : lHildebrandt et al.1l2007l : iGold et al.ll2008l ). Given that 
the DL98b models are now a decade old, and the recent surge 
in interest in anomalous emission, it is timely to revisit the 
theory of spinning dust emission, including the approxima- 
tions made in DL98b. This is the purpose of this paper. 

As in DL98b, we concentrate on the rotation rate of the 
grains; the size distribution has been reconsidered by other 
authors, and the grain dipole moment distribution should 
be regarded as a model parameter since one cannot com- 
pute it from first principles. We first review and generalize 
DL98b rotational excitation and damping rates. We mod- 
ify the rotational excitation and damping rates by collisions 
with neutral species, such that it respects detailed balance 
in the case where the evaporation temperature is equal to 
the gas temperature. We include the electric dipole potential 
when evaluating the effect of collisions with ions. Full hyper- 
bolic trajectories and rotating grains are used when comput- 
ing the effect of plasma drag. We correct the infrared emis- 
sion damping rate which was underestimated for a given in- 
frared spectrum. Finally, we use these excitation and damp- 
ing rates to calculate the grain rotational distribution func- 
tion by solving the Fokker Planck equation. Updated grain 
optical properties and size distribution are used throughout 
this analysis. An Interactive Data Language (IDL) code im- 
plementing the formulas in this paper, SpDust, is available 
on the wel0, and will hopefully allow for a more thorough 
exploration of the parameter space, as well as model fitting 
to observations. 

The paper is organized as follows. In Section [2] we re- 
mind the reader of the electric dipole radiation formula and 
give the resulting expected emissivity. In Section [3] we dis- 
cuss the size distribution and dipole moments, along with 
other grain properties. We then turn to the main thrust of 
this study, which is the computation of the angular velocity 
distribution function. The theoretical formalism is exposed 
in Section |4j which presents the Fokker-Planck equation. 
Sections [SHI] discuss the various rotational damping and ex- 
citation processes : collisions with ions and neutral species, 
plasma drag, infrared emission, photoelectric emission, and 
random H2 formation. The reader interested primarily in 
the predicted emission may wish to proceed directly to Sec- 
tion 1101 where we present the resulting emissivity and the 
effect of various parameters and environment conditions. 
Our conclusions are given in Section [11] Appendix A ex- 
poses the techniques used to numerically evaluate integrals 
of rapidly oscillating functions involved in the plasma drag 
calculation. Appendix B presents an alternate, quantum me- 
chanical derivation of the rotational damping rate through 
infrared emission. 



2 ELECTRIC DIPOLE RADIATION 

The power radiated by a dust grain spinning with an angular 
velocity uj, of electric dipole moment /i, with component /ix 
perpendicular to aj, is 



^ http:/ /www. tapir. caltech.edu/~yacine/spdust/spdust. html 



This power is emitted at the frequency u — uj/2n. 

To get the emissivity of electric dipole radiation per H 
atom, in ergs~^ sr~^ (Hatom)"^, one needs several ingredi- 
ents: 

• The grain size distribution function: n^^dugr/da, which 
gives the number of dust grains per unit size per H atom. 

• The electric dipole moments as a function of grain size 

a: 

• The angular velocity distribution function, fa {i^) , which 
depends upon the grain radius and environmental condi- 
tions. It depends on the angular velocity modulus only in 
a perfectly isotropic environment, with no strong electro- 
magnetic fields forcing the dipole moments to align in some 
particular direction. 

One then readily gets the emissivity of spinning dust grains 
per H atom: 

= — / da — 4ttuj fa{uj) 2n — , (2) 

riH 47r Ja„i„ "■H da 3 cP 

where ui = 2-kv. 



3 DUST GRAINS PROPERTIES 
3.1 Grain shapes 

The grains are characterized by their volume-equivalent ra- 
dius a, such that the grain volume is 47ra^/3. The radius a 
is in fact a measure of the number of C atoms in the grain, 
which we assume to be 

= 1^ ^ 468 al, (3) 
3mc 

where pc = 2.24g cm~^ is the density of ideal graphite and 

a_7 = a/(10~^ cm). 

We foUow iDraine fc Lil (|200lh . hereafter DLOl, for the 
number A'h of H-atoms in the grains (see their Eq. 8). Fol- 
lowing DL98b, we account for the fact that the smallest 
grains may be sheetlik^l, as expected for polycyclic aromatic 
hydrocarbons (PAHs). We assume that this is the case for 
a < a2 = 6 A (this corresponds to A'^c ~ 100 carbon atoms, 
the size of a large PAH). We model them as disks of thick- 
ness d — 3.35 A, the interlayer separation in graphite. In 
many cases, these grains will be rotat ing primarily around 
the axis of largest moment of inertia (jPurce lil ll979l ), which 
is perpendicular to the plane of the grain. When comput- 
ing various rates, we will usually assume a spherical geom- 
etry, with a "surface-equivalent" radius as or a "cyhndrical 
excitation-equivalent" radius Ucx, defined as : 

47ra^ = ^ dS and 47raL = ^ ^ P^dS', (4) 

where p = rsinO is the distance to the axis of symmetry 
and dS is the surface area element. 

Although the assumption of cylindrical grains for a < a2 
is not critical, it does have an effect on the spectrum, which 
is shown in Fig. [1] 

^ DL98b allow for a possible population of linear grains, although 
they do not actually use them. 
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Cylindrical if a < a, = 6.10 cm 
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Figure 1. Effect of relaxing the assumption of cylindrical grains 
on the spectrum, for a fiducial Cold Neutral Medium environment 
(CNM, defined in Eq. I|162|l ). At equal radius, spherical grains 
have a smaller moment of inertia than the cylindrical ones, which 
are rotating primarily about their axis of largest moment of iner- 
tia. They will thus radiate at slightly higher frequencies. For the 
CNM, we find an increase of peak frequency A^'peak/!^peak ~ 6%. 
The high-frequency tail of the spectrum is due to the smallest 
dipole moments of the assumed Gaussian distribution for the in- 
trinsic dipole moments (see Section l3.3l and Fig llOH . For a spher- 
ical distribution of dipole moments, there are fewer grains with 
a low dipole than for a planar distribution. This explains the de- 
crease in power at high frequencies. For the CNM, this results in 
a decrease of total emitted power Ajtot/itot ~ —16%. 



3.2 Size distribution 

Following IWeingartner fc Draind (j2001al l. hereafter WDOla, 
we adopt the following size distribution for carbonaceous 
dust, for grain radii amin = 3.5 A< a < amax = lOOA: 



1 drigr 
TiH da 



where 



F{a;f3, at) 



Dia) + ^UL 
a \ at 



1, 



-l{a~at)/a^]^ 



F(a-p,at) 

< a < at 
a > at 



1 -I- Pa/at, 
(1 - Pa/ at) 



13 <0 



(5) 



(6) 



The function D{a) describes truncated (at 3.5 A) lognormal 
grain populations, 



D{a) = ^ — exp 



ln(a/ao,i 



(7) 



with the normalizations Bi defined to place a total number 
bc,i of carbon atoms per hydrogen nucleus in the i"' log- 
normal population. Here 6c, i = 0.75&C, be, 2 ~ 0.25bc, &c 
being the total carbon abundance per hydrogen nucleus in 
the lognormal populations, ao,i — 3.5 A, ao,2 = 30 A, and 
a = 0.4. This size distribution has a total of six adjustable 
parameters {be, C, a^, a^, a, 13). For a given choice of be, the 
other five parameters can be found in WDOlb, Table 1. 

We consider only carbonaceous grains in this work. The 
abundance of small silicate grains in the diffuse phases is in- 



deed limited by the absence of the 10 /im band in emission, 
as disc ussed in WDOlb. Note, however, that iLi fc Draind 
( :2001al ) found that as much as ~ 10% of the interstellar 
silicate mass could be in the form of ultrasmall particles 
(a < I5A) without violating any existing observational con- 
straints. While consistent with observations, our assumption 
is thus not required by them. 



3.3 Dipole moments 

Although only the component of the dipole moment perpen- 
dicular to uj is of importance for the electric dipole emission, 
the total dipole moment will be needed in coming calcu- 
lations. In case of spherical grains, we assume the dipole 
moment and rotation axis are randomly oriented. For cylin- 
drical grains, the dipole moment is mainly oriented in the 
plane of the grain, perpendicular to the rotation axis. 

The dipole moments have two components. First, an 
intrinsic part /Xi, which results from the addition of dipole 
moments from individual molecular bonds. For a given grain 
radius, we assume a multivariate Gaussian distribution, with 
variance proportional to the number of atoms in the grain, 

iVat =Nc+Nh- 



P{pi) oc 



with 



spherical grains 
disklike grains 



(8) 



(9) 



These are appropriate assumptions if the dipole moments of 
bonds add in a random-walk fashion, although we caution 
that this need not be the case. Counterexamples could in- 
clude dipole moments dominated by a single feature, e.g. a 
PAH with a single OH group. The formula given above is 
in that case intended to give an estimate of the total dipole 
moment, with the value of /? tuned to reproduce approxi- 
mately observed dipole moments for laboratory molecules 
(see e.g. DL98b, Table 3). 

The rms dipole moment per atom, /3, is poorly known. 
Following DL98b, we will take /3 = 0.38 debye as a fiducial 
value, corresponding to 



(M?)la = 



10-'' cm 



= (9.3 debye) 



(10) 



In addition, for grains with charge Ze, a displacement d 
between charge centroid and center of mass (e.g. due to 
asymmetric grain shape or isotopic substitution) may add 
another, uncorrelated component. We assume that the dis- 
placement is proportional to the excitation equivalent ra- 
dius: d = ea^x, where e = 0.01 (DL98b). In most cases this is 
negligible compared to the intrinsic component, so we model 
it as a single value for the sake of simplicity. The total dipole 
moment is thus given by 

^'^ = fii + {eZqcaax)^ , (11) 
where qe is the elementary charge. 



3.4 Grain charge 

The rotational damping and excitation rates will be depen- 
dent on the grain charge. DL98b showed that the character- 
istic timescale for changes in charge is much shorter than the 
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characteristic rotational damping time. We will therefore av- 
erage the damping and excitation rates over grain charges, 
as well as the electric dipole moment when computing the 
power radiated. We therefore need the charge distribution 
functiorffl of the grains as a function of their radius and en- 
vironmental conditions, fa{Z). 

There are three main processes contributing to grain 
charging: coUisional charging by electrons and ions, which 
rates we denote Je{Z,a) and Ji{Z,a) respectively, and pho- 
toelectric emission of electrons caused by the impinging ra- 
diation, which rate is Jpe{Z,a). For every grain radius, the 
steady state charge distribution function is obtained by solv- 
ing recursively the following equations: 

[MZ, a) + Jp4Z, a)] fa{Z) ^MZ + 1, a)fa{Z + 1). (12) 

We use the equations of iDraine fc SutinI ll 19871) for coUi- 
sional processes, updated with the Weingartner fc Draind 
l|2001bl ) electron sticking coefRcients, for Ji and Je. The 
photoelectric emission rate is computed according to 
WDOlb. The radiation field is taken to be a multi- 
ple X of the ave rage interstellar radiation field misrf, 
as e stimate d by Mezger, Mathis fc Panagial (|l982l ) and 
iMathis. Mezger fc Panagial (|l983l ). 



4 THE FOKKER-PLANCK EQUATION 

4.1 Form of the equation in spherical polar 
coordinates 

The stationary angular velocity distribution function fa{^) 
(such that /a(<j;)d"^a; is the probability of the grain's angular 
velocity being u: within d^tj) is determined from the station- 
ary Fokker-Planck equation. We differ here from DL98b who 
assumed the distribution was Maxwellian and calculated its 
approximate rms grain rotation rate < lo^ >^^^ . The Fokker- 
Planck equation is valid in the limit of continuous torques, 
i.e. if every interaction changing the rotation rate of the 
grain does so by a small amount Suj <^ uj. This is, therefore, 
accurate for the largest grains, which have large moments 
of inertia. But it fails to describe precisely the smallest ones 
(a < 7 A), for which DL98b showed that impulsive torques 
are important (see their Section 7 and Fig. 7). However, 
we believe that the actual distribution function would dif- 
fer from the one we calculate only at very high frequencies, 
where the dust emissivity is dominated by the vibrational 
emission. Indeed, the occasional impulsive torques on the 
grains enhance the distribution function for high values of 
the rotation rate, where the solution of the Fokker-Planck 
equation predicts an exponential cutoff, as we shall see later. 
The peak of the distribution will not be affected significantly, 
as the variations of the rotation rate of a grain within the 
peak are not impulsive. 

The stationary Fokker-Planck equation is given by 

^ [d^(.)/.(c.)] + = 0. (13) 

The coefficients are defined as: 



^ We use the same notation for different distribution functions. 
The context and their argument should make their meaning un- 
ambiguous. 



D\uj) = - lim ^— ^ and E'\uj) = lim ^ ; ^.(14) 



st^o St 



We assume that the medium is isotropic, and there are no 
physical processes that allow for a preferred direction, such 
as a magnetic field. As a consequence, the rotational distri- 
bution function only depends upon the magnitude uj of u>. 
Moreover, in a local orthonormal frame {eu,,e0,e^), where 
Lo, 9 and (j) are the usual spherical polar coordinates defining 
u:, the excitation coefficient take up the following form : 



accounts for fiuctuations along lJ, and 



(15) 



(16) 



accounts for fiuctuations perpendicular to Ct). The compo- 
nents in the coordinate basis are thus: 



E--^E^^iLo), e''^^^, e^^ 



E^{lo) 



(17) 



Moreover, we assume there are no systematic torques, so the 
damping coefficient is directed along lj and we have 



(18) 



In the spherical polar coordinate basis, the Fokker-Planck 
equation then becomes: 

[c.^i?ii(^)A(^)] 



[uj E^{uj)Uuj)] 



0. 



(19) 



Integrating once, we get the following first order differential 
equation: 



where 

1 , ,1 d£;|| 

LO 2 dcj 



(20) 



(21) 



Note that D is simply equal to D if the fluctuations are 
isotropic and independent of ui. 

The coefficients D, E±, and therefore D from var- 
ious independent rotational damping and excitation pro- 
cesses are additive. 

A given process is said to respect detailed balance, 
when, if that process were the only one taking place, the 
grain would rotate thermally, i.e. fa{i-o) oc exp{—Iuj^/2kT). 
As one can see from the Fokker-Planck equation, this implies 
that this process must satisfy : 



D 



2kf' 



En 



(22) 



Excitation rates are often easier to calculate than damping 
rates, since they are positive definite and do not rely on near- 
cancellation of processes that increase versus decreasing lu. 
Thus in some cases, we will make use of detailed balance (i.e. 
the fiuctuation-dissipation theorem), to obtain the damping 
rate, knowing the excitation rate. 
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4.2 Normalized damping and excitation 
coefficients 

We will see in the next section that for collisions with neu- 
tral H atoms, at a temperature T, for a spherical dust grain 
at the same temperature T, the damping and parallel exci- 
tation coefBcients have the following form: 

2kT 



Dh ^ — and E\\_„ = E\m = 



TH 



where 



TinmH 



2fcT 
Trmn 



1/2 



3/ 



(23) 



(24) 



is the characteristic rotational damping timescale for col- 
lisions with neutral H atoms. Note that they respect the 
detailed balance condition. 

We normalize the damping and excitation coefBcients 
of each process to those of collisions with H atoms. Taking 
DL98b notation, we define, for each process X : 



Fx{uj) 



(25) 



Gx(a;)= ^£||,x(^) (26) 

A special case is made of the rotational damping through 
electric dipole radiation (subscript ed), because of its specific 
OJ^ dependence: 



dt v^^"' 



so 

Auj 
dt 



2 ti\uo''^ 



3 /c3 

Here we define, following DL98b: 



3kT Ted 



Ted 



2fcr/i2 



(27) 



(28) 



(29) 



Using Eqs. ([251, ^ and ([28} in Eq. (O, the final equation 
for the distribution function is 



dfa 



+ 



2, ,3 1 



1 LU 



where 



Icq F th 1 

kTG ^ ^3G {kT)2 



(30) 



X and Gee^Gx. (31) 

X X 

One can see that the conditions to get a thermal, Maxwellian 
distribution fa((^) oc exp{—Iuj^ /2kT) are: 

F = G = constant and > 0. 

Ted 

Otherwise, the general solution to this equation is 

luj' F{u') 



(32) 



fa{i^) OC exp 



duj' 



kT G{uj') 



+ 



3redG(aj') (fcT)2 
If all Fx's and Gx's are constant, this has a simple form 



(33) 



fa{i^)(x exp 



Fiij_ _ TH 1 fi'^^y 

'G2fcT ~ ^3G V2fcT/ 



(34) 



Note that the damping through electric dipole radiation 
causes the distribution to be non-Maxwellian. 

In the general case, some Fx's and Gx's may depend 
upon Lo and one has to compute numerically the resulting 
distribution function, using Eq. p3p . 

We now turn to the calculation of the various damping 
and excitation coefficients, due to collisions, plasma drag, 
infrared emission, photoelectron emission, and random H2 
formation. In the following microphysics sections that form 
the heart of the paper, we compute excitation and damping 
coefficients as a function of grain radius and environmental 
conditions. We evaluate them numerically for a fiducial Cold 
Neutral Medium (CNM) environment, defined explicitly in 

Eq. (Hsg. 



5 COLLISIONAL DAMPING AND 
EXCITATION 

In this section we correct the results of DL98b, Appendix B, 
which did not take into account the fact that not all neutrals 
escape the grain surface when computing the damping rate. 
The microphysics of collisions is complex and beyond the 
scope of this study (for a discussion of the physics and chem- 
istry of PA Hs and their re lation with the interstellar gas see 
for example lOmontI (|l986l )'). We therefore use the following 
simplifying assumptions: 

• The grain is in a stationary state: the rate at which 
species collide with it is equal to the rate at which they 
leave its surface. 

• We assume that all species (neutrals and ions) colliding 
with the grain stick and that they depart the grain as neu- 
trals. In extremely dense environments, the colliding species 
may bouce off the grain surface instead of sticking. This case 
is discussed at the end of Section 15.1.41 

• Even if the impacting species may not collide equiprob- 
ably everywhere on the grain's surface (e.g. if the grain is 
non spherical or if it has a dipole moment), we assume they 
somehow get re-distributed on the grain surface and leave it 
equiprobably from any point. 

• We assume, as in DL98b, that neutrals leave the grain 
surface with a thermal velocity distribution in the grain's 
frame, with a temperature T^v of the order of the infrared 
emission characteristic temperature. Unlike DL98b, we esti- 
mate Tev as a function of grain radius and ambient radiation 
field (see section r5.1.4|) . 

Using those assumptions, one can compute the rate of col- 
lisional damping and excitation. We will perform the calcu- 
lations for a spherical grain in the general case. To find the 
relevant equivalent radius to use for a cylindrical grain, we 
will carry out the explicit calculation in the case of colli- 
sions of a neutral grain with neutral H atoms. Note that as 
pointed in DL98b, the rotational excitation in case of col- 
lisions has two origins: the random excitation by incoming 
particles (superscript '^*"'), as well as the random excitation 
by "evaporating" neutrals (superscript ('=")). 
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5.1 General considerations: spherical grain 



get: 



We use the usual spherical polar coordinates around the 
spherical grain, taking the rotation axis as a reference. The 
local phase-space density at the grain surface is: 



with the local velocity 



: uj X r = au) sin 9e 



0. 



(35) 



(36) 



The normalization constant K(6) is found by imposing that, 
locally, the flux of evaporating (and escaping) particles is 
equal to the flux of colliding particles. Except for the case 
of ions interacting with the electric dipole of the grain, the 
flux of colliding particles will be homogenous on the grain 
surface. If it is not the case, we approximate the local flux 
by the total rate of collisions dNcoii/dt divided by the grain 



1 dN,M 
Aiza^ At 



K / Vr CXp 



V 



-Pcscd V, 



(37) 



where Pose ~ 1 for velocities at the grain surface leading to 
escape, and otherwise. 

All particles evaporating from the grain are neutrals. 
They interact with the grain through the induced dipole 
potential (we neglect the dipole-induced dipole interaction 
with the dipole moment of the grain): 



(38) 



where a is the polarizability of the escaping neutral and 
qe is the elementary charge. The polarizability of hydro- 
gen is a standard result in nonrelativistic quantum me- 
0.67 where ao is the Bohr radius 



chanics and is ^af, 



l|Landau fc Lifshit j|l965l). We also ta ke a = 0.20A^ for he- 



liurrQ (jThomas fc Humbcrt sonlll972l 'l. and a = 1.54A^ for 



carbon (Miller & Kelly 197^ which is important since C 
is often the dominant ion if the hydrogen is se lf-shielded. Fo r 
molecular hydrogen H2, we take a = 0. 79A^ (|Marlowlll965h . 

5.1.1 Computation of Pesc 

The radial coordinate of the escaping neutral is the solution 
of the following equation: 



2E 



r + Vcff(r) = r + -irVu jVa = 



(39) 



where v\\ is the modulus of the tangential velocity at the 
grain surface and 



2 _ ^s^eO 



The effective potential has a maximum at the radius 



(40) 



(41) 



To escape, a neutral needs to have either a > or 

2E/m > Vcs{ra). These two conditions can be combined to 



We assume that all the helium is neutral and riHo/'i-H = 1/12 



Pose = 1 if 



Vr > Va or 



< Vr < Va and lly > \/2Va(Va — Vr) 

(42) 

where Vr is the radial velocity at the grain surface. 
5.1.2 Computation of K (6) 

Following DL98b, we define = mVa/2kTev, which de- 
scribes whether the typical evaporating atom has enough 
energy to overcome the induced dipole attraction to the 
grain (ce < 1) or not {te > 1). We also define the ratio of 
rotational velocity to thermal velocity at the grain surface, 
which is small compared to unity: 

1/2 



O, = auj 



2kTe', 



< 1. 



(43) 



In terms of those dimensionless quantities, we can find the 
normalization constant K. The right-hand side of Eq. (|37|l 
can be expanded using the substitution 



{Vr,V0,V^) = 

to yield 

1 diVcoii 
47ra2 At 



{ur^u COS ijj^ u sin -i/;) 



(44) 



K 



^k'J. (=■] 



-(Hsin^)^ 



2UrAUr e 



2uAue " Io(2uilsm9) 



-Ur) 



where 



1 

2^ 



X sin ^ J I 

e dip 







1 + -X^ + ... 

4 



is a modified Bessel function of the first kind 
Expanding to second order in fl, we get 
-2 „ ,2 



K = 



V 



1 rfiVco!, 



(45) 



(46) 



(47) 



\ m J TT e~^i + 0ree erf (ee) 47ra2 dt 
up to corrections of order 0{Q,'^). 

5.1.3 Damping and excitation rates 

Each escaping neutral particle takes away an angular mo- 
mentum 



L = ma{vee^ — v^ee). 



(48) 



As Pesc is an even function of vg , the average of vg vanishes. 
The loss of angular momentum along the z-direction per unit 
time per unit area is given by 

-ma sin ^ K I UrU^ exp 



AL, 



At AS 



mlv — Vo ) 



XPeac AVr Avg Av^. 



(49) 



Here we differ from DL98b as we take into account the fact 
that not all particles escape from the grain. Expanding in SI 
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and using the expression for K we get, up to corrections of 
order 0{D.'^): 



dtdS 



1 . 2, 

msm y — 

47r e- 



dN„ 



+ ypKee erf dt 



(50) 



Integrating over the whole grain surface, we find the damp- 
ing rate 



IdL^ 
7~dr 



2ma^ diVc< 



+ \/Tree erf te 3/ dt 



(51) 



A similar calculation leads to the excitation rate through 
evaporating particles: 

_ e"'= + 2el 2ma^ dA^coii 

e-^e +y;Feeerfee~372 d^ 



up to terms quadratic in Q. 

This implies the remarkable relation 



G 



(ev) 



2T' 



■Fc< 



(52) 



(53) 



Physically, this occurs because if Tev = T then the collisions 
with neutrals satisfy detailed balance, Eq. (|22p . The factor 
of 2 arises since in this case there is an equal contribution 
to the excitation from incoming and evaporating particles. 

We derive a stronger damping rate due to evaporating 
atoms than DL98b: for ^ 1 this results in no change, but 
for Ee ^ 1 we find much stronger damping. The physical ori- 
gin of this is that atoms that evaporate with prograde veloc- 
ities relative to the local grain surface {v^ > vori,) typically 
have more angular momentum than atoms that evaporate 
with retrograde velocities. Therefore the centrifugal poten- 
tial helps them to escape the grain. DL98b neglected this 
effect, but for ee S> 1 it is dominant. 

The excitation rate through incoming particles will be 
calculated for each case. 



5.1.4 Evaporation temperature Tev 

DL98b assume that the evaporating temperature is a con- 
stant, independent of grain size. This accurately describes 
the largest grains, for which the temperature may be approx- 
imated as a constant, obtained from equating the absorbed 
and emitted energy (DL98b): 



he 
T 



{Q),u, 



87r/ic(3oA^r(a + 4)C(a-|-4) 



(54) 



where in the infrared the grain absorption efficiency is as- 
sumed to be a power-law 



= Qo 



Ao = 



1^0 



(55) 



with typically a — 2 and {Q),Ui, = J di/Q^u^. Note that we 
have Tc tx x^''^ and a weak dependence on grain radius as 
the absorption efficiencies cancel out. 

However, the smallest grains undergo sudden thermal 
spikes after each photon absorption, followed by long inter- 
vals during wich the grain drops to its vibrational ground 



1000 r 




a (cm) 



Figure 2. Evaporation temperature T^v as a function of grain 
radius a, for various values of the ambient radiation field u, pa- 
rameterized by X = f/flSRF- The curves join at small radii, 
for which the grains undergo temperature spikes. The kink at 
a = 50A results from the DLOl prescription for PAH-graphite 
optical properties. 



state. The neutrals or ions that have stuck to the grain after 
a collision cannot be thermally ejected from a grain in the 
ground state so we assume ejection during thermal spikes. 
A simple assumption is that in this case ejection occurs af- 
ter a photon absorption and thermalization of the photon's 
energy. We take 



/ QvUi^dv 



(56) 



as the typical energy of an absorbed photon. Typically, E.^ ~ 
5 eV. We then calculate the corresponding grain tempera- 
ture following DLOl: we solve for Tq such that E{Tq) = E.^, 
where 



tlLOi 



ixp{riujj/kT) - 1 



(57) 



is the expectation value of the energy of the grain, and the 
sum runs over its Nm vibrational degrees of freedom. We 
take Tev = max(rc,T,) as the evaporation temperature. 
The result is shown in Fig. [2] One can see that we obtain 
much higher evaporation temperatures than the ones used 
by DL98lJ!|. The effect may be significant on the final spec- 
trum, as can be seen from Fig. [3] 

High density, low radiation field case 

The previous treatment is valid only if the rate of photon 
absorption is high enough to eject all stuck species before 
all available sites on the grain are occupied. We approximate 
the number of available sites on the grain by the number of 
superficial C-atoms : 



^ The mechanism we describe for atomi c ejection from gra ins is 
called photo-thermo-dissociation (PDT). lRouan et al.l il99t i also 
mention another possible mechanism, photo-dissociation (PD), 
which is an atomic ejection following the direct interaction of a 
UV photon with a given C-H bond. PD may lead to even higher 
ejection temperatures, of order 10 000 K 
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V (GHz) 



Figure 3. Effect of the evaporation temperature model on the 
spinning dust spectrum for the Cold Neutral Medium (CNM, Eq. 
I|162|l ) . Our prescription results in a much higher evaporation tem- 
perature for the smallest grains, compared to DL98b, who assume 
a constant T^v = 20 K for all grain sizes. This leads to a decreased 
damping rate (see discussion at the end of section fS. 1.31 1 and an in- 
creased excitation rate through collisions, and therefore increases 
the peak frequency of the spectrum. 



iVsitc 



Nc{a) 
MiVc(a) 



for cylindrical grains 
for spherical grains 



(58) 



where Ncifl) was defined in Eq. ([3]) and d — 3.35A is the 
interlayer separation in graphite. The ratio of collision rate 
to photon absorption rate is given by : 



'coll/abs 



TTTTlH 



-du c 



' O.lx- 



riH 



30cm 



In most environments, /?coii/abs «^ A'sitcs so there is no ac- 
cumulation of stuck species. In very dense and dark clouds 
however, the rate of collisions may become so high com- 
pared to the rate of photon absorption that all the sites 
are occupied. In that case, the assumption that incoming 
species stick to the grain is no more valid. They will instead 
bounce off the irregular grain surface. From the fluctuation- 
dissipation theorem, one expects that, for collisions with 
neutral species, F„ = Gn- Thus, we set the effective evap- 
oration temperature equal to the gas temperature in that 
case (see Eq. H74p and discussion below): 



Tev =T if i?e oll/abs 



(60) 



The actual transition from sticking to elastic collisions 
should of course be smooth, unlike the discontinuous step 
we assume here. Our treatment should approximately reflect 
the physics of collisions except near the transition regime 

-Rcoll/abs ~ A'sitcs. 



fin{v) 



3/2 



-mnf /2feT 



(61) 



\2TTkTJ 

from which one can easily get the excitation rate through 
incoming particles : 

r 2 



dALj 
dtdS 



= / v„{miipv^) fi„{v)d V 



(62) 



where ii„ is the component of the velocity normal to the 
grain surface. This evaluates to 



2kT 



3/2 



Integrating over the grain surface, we get 



dAL 

dt 



2 {»") 



= kTriurriH 



2kT 



1/2 A 4 



(63) 



(64) 



TTTriH J 3 

where acx was deflned in equation (j4|. For a spherical grain, 
Ucx ~ a. For a disklike grain of thickness d and radius b, 
spinning around its axis of symmetry, we have 

-,1/4 

(65) 



lb' (2d + 6) 



We can write the excitation rate by incoming H atoms as 
kT 



j-,(in) 



Ith 



(66) 



where th was deflned in equation (|24|> . 

The case of evaporating particles is very similar. As- 
suming the grain surface is at the same temperature T as 
the gas, the phase-space density of evaporating particles is 



.Uv{v) 



( triH \3/2 

V27rfcr/ 



exp 



2kT 



(67) 



In that case P^sc ~ 1 for all outgoing particles. The same 
calculation therefore leads to 

kT 



-^11, H - 



/th 



(68) 



up to terms of order 0{Q^), which comes from the fact that 
we did not take into account the slight change of ui after 
the particle has collided (we assumed the same to for the in- 
coming and the outgoing particle). Detailed balance ensures 
that 



TH 



(69) 



Therefore, for non spherical grains, we will compute collision 
rates assuming a spherical geometry with radius acx ■ We just 
showed that this is an exact result for collisions with neu- 
tral H atoms. The collision rates are indeed proportional to 
the area of the grain, but the angular momentum gained 
depends on (p^), so acx will approximately reflect both de- 
pendencies. 



5.2 Collision with neutral H atoms: neutral grain, 
general grain shape 

We assume that the grain is neutral, and has no dipole mo- 
ment, so there is no interaction whatsoever between the 
grain and the neutral H atoms (purely geometric cross- 
section). The phase-space density of incoming H atoms at 
the grain surface is simply 



5.3 Collisions with neutral atoms: charged grains 

In that case, the incoming neutrals interact with the same 
potential as the outgoing particles: 

U{r)^-\a^. (70) 
We use the same notation as DL98b and deflne 
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>— and 6oW = aW— , 



(71) 



where Va was defined in Eq. (|40p . 

We recall, from DL98b, that a trajectory with impact pa- 
rameter b and velocity at infinity v leads to a collision if 



b 4: 6max(u) 



foo(«) if V !^Va 

a^l + vl/v'^ if n ^ 



We compute the collision rate 



diVc, 



df 



( ran 



i/2 



(72) 



-mv^ /2kT 



nn2na 



(2kT\ 
\ Tim J 



1/2 



\2TikT) 
l^e""" + \/7re„ erf e„ 



(73) 



We can now get the normalized damping and excitation 
rates for collisions with neutrals: 



Fn 



nn I m„ e ^" + y/iren erf e„ ( 2\ 



G 



(ev) 



riH V e^'^e -)- yTTEe erf e, 



2T 



Fn, and 



G, 



(in) 



2ni^ Y mn V "/ 



(74) 



where the result for Gil"'' is identical to that of DL98b. 
Note that when T = Tev, c'n"^ = gI™' = F„/2 so the 
principle of detailed balance holds. Moreover, in the case 
of a neutral grain, if the only rotational excitation and 
damping process were collisions with neutral species, then 
the rotational distribution function would be a Maxwellian. 
In that case the rotational temperature would be given by 
Trot =Gn/FnXT = ^{T + n^), the arithmetic mean of the 
gas and evapo ration temperatures, as was already shown by 
IPurcelll (|l979l ). 

This is the contribution of an individual neutral, for 
a given grain charge. To get the total contribution, one 
must average over all grain charges (DL98b showed that 
the charging timescale is much shorter than the collision 
timescale) and sum over all neutral species, which we take 
to be atomic and molecular hydrogen, and heliurrQ (with 
nuc/riH = 1/12). 



5.4 Collisions with ions: charged grains 

The ion interacts with the grain through the Coulomb, 
electric dipole, and "image charge" potentials. The latter 
dominates over the Coulomb potential only in the imme- 
diate vicinity of the grain surface, so we will neglect it for 
charged grains. Properly accounting for it would result in 
a slight increase in both damping and excitation rates as 
this is an attractive potential. The general solution for this 
problem, with a rotating electric dipole moment, is still not 
analytical. Thus, for simplicity, we will only consider the 
case where the electric dipole moment can be considered 
as non-rotating, i.e. when the timescale of the collision is 



° Collisions with neutral helium have little effect on the spectrum: 
the helium contribution dominates Fn and Gn only in the case 
where the medium is strongly ionized, i.e. when the dominant ro- 
tational excitation and damping mechanisms are rather collisions 
with ions or plasma drag. We include them for completeness. 



short compared to the rotation period of the grain. This is 
justified as, when the ion reaches the vicinity of the grain 
surface, the ratio of collision timescale to rotation timescale 
is approximately Lva/v ~ ^rrii/mgr <^ 1. We will assume 
that the grain is spherical, so that the electric dipole 
moment is randomly oriented relative to the rotation axis 
(for cylindrical grains this is not the case but we will assume 
so for simplicity). Note that when the grain rotates rapidly, 
the component of the dipole moment perpendicular to the 
rotation axis averages out, but not the parallel component. 
Although this problem will be different in nature as this 
alignment creates anisotropic excitation by collisions, the 
magnitude of the non rotating part of the dipole moment 
will remain of the same order (a factor smaller only), 

so our approximation should give a decent idea of what the 
effect of the dipole moment is on the trajectory. 

We assume a spherical geometry with radius a^x- Tak- 
ing fj, as the polar axis for spherical polar coordinates, the 
interaction potential of the ion in the Coulomb and dipole 
field of the grain is given by 



V{r 6) = ^g^'*^^ + Ziqencose 



(75) 



r 

The motion in this potential has two obvious constants: the 
energy E and the angular momentum along the jz-axis (along 
fi), Lz. For this special potential, however, there exists a 
third constant of the motion. The torque L exerted on the 
ion comes entirely from the second term in the potential and 
is 



L = -r X \7V{r) = — e,^ 



(76) 



Since the azimuthal component of angular momentum is L ■ 
etf, = mr^6, we can then determine the overall rate of change 
of the angular momentum, 

^(L^) = 2miZiqcfi sin 9 6 = —2miZiq^,^^cos9. (77) 
Therefore we find the constant of the motiorQ 



A = L + 2miZiqcficos( 



(78) 



Its value can be determined by taking the incoming trajec- 
tory at infinity with approach angle 9oa, 



A- ~ (rriibv) + 2miZiqc^cos9oo 
The energy of the trajectory is 
1 2 

-mj-r + T4ff(r) = E 



(79) 



(80) 



where Vefi(r) is the sum of the potential V(r,9) and the 
tangential kinetic energy /2mir'^: 



■tr r \ - ZgZiql mt{bv)'^ + 2Ziqcfi cost 
Vcs(r) = 1 



2r2 



(81) 



It is easier to work with the following dimensionless param- 
eters 



ZgZjq^ 

kT 



(82) 



^ This may also be derived by the Hamilton-Jacobi method in 
spherical polar coordinates. 
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Their physical meaning is as follows : \tp\ <^ 1 when the 
thermal energy of the ion dominates over the electrostatic 
interaction energy, and [tpl ^ 1 when the electrostatic in- 
teraction dominates. The sign determines whether the inter- 
action is attractive {tp < 0) or repulsive (i/j > 0). /i is the 
equivalent quantity for the dipole interaction. Note that we 
consider only positively charged ions so jl > 0. We also work 
with the dimensionless variables 



U — \ I ——;V. 

2kT 



The effective potential can now be written 



Kff(r) = fcT 



^ }- [u c + jl cos f 



(83) 



(84) 



A study of this potential leads to the following condition for 
collision : 



< Xmax(c, u) = /i [u -UC - 1p) 



(85) 



Note that if Xmax < —1, then there is never collision, for 
any angle. If Xmax > 1, then all angles lead to a collision. 
We define 



X{c,u) = max{-l,min [1, Xmax(c, u)]} 
Now, we can compute the collision rate 
dA^coii 



(86) 



At 



m f 2nv^dv 2TTbdb ( 

J \2-KkT 

xe-™'"'/2'=^(X + l) 

l2TxkT 2 

2ni\ a^^ 

V rrii 



3/2 



X / 2M''e"'' Au 2c dc 



X + 1 



We can also get the excitation rate by incoming ions 



(87) 



dt 



{mhvY 



2-Kv-'Av 2TvbAb 



\2nkT 



3/2 



xe-'"'^'/2'=^(X + l) 



2nimi a^^ ( 2-RkT 



3/2 



Stt V rrii 
X / u°e~''~ Au Ac' Ac 



These integrals can be evaluated explicitely and one then 
gets, for the charged grains Zg 




Gf"'(Z,/0) 



o 1/ — 5'2(^,m), 
2nH V "^H 



(89) 



where we have defined g\{xp,ji) — 

1-tp V < 

e"''' sinh fi/fl tp > 

2/i 



(90) 



and g2{tp,il) = 

e'"^ smhfl/fl ip > 

^-Ki(M-^)' + KM-'^)' 



2/1 



, M> 1^1- (91) 



Note that these functions coincide with the functions 
gi{ip) , g2{tp) defined in DL98b for /i = 0. We also defined 
6; = Zgq^ai/2a'^kTev (here a; is the polarizability of species 
i after it neutralizes on the grain surface, e.g. when consid- 
ering collisions with the C"'" ion we take the polarizability of 
the neutral C atom). Note that even when Tei, — T, Fi d 
as the incoming and outgoing particles are in different ion- 
ization states; detailed balance does not apply since realistic 
ISM phases are not in Saha equilibrium. Numerically, one 
has (with T2 = T/IOOK) 

(92) 
(93) 



-tp ^ 170 Zg aZl Ti^ 



30 



2\1/2| 



'cm -1/2 rr-1 



9.3 debye 

From these values, one can see that in general the effect of 
the dipole moment cannot be neglected a priori, as jl is not 
small compared to unity. However, in general fi < \tp\. This 
implies that, for negatively charged grains, the dipole mo- 
ment has little or no effect on the excitation and damping 
rate. For positively charged grains, the damping and excita- 
tion rate are both increased by the huge factor sinh /t//i, but 
still remain extremely small due to the Coulomb repulsion, 
which shows in the factor e"'''. 

We therefore conclude that DL98b approximation of ne- 
glecting the effect of the electric dipole moment on the tra- 
jectory of ions, is essentially valid in the case of collisions 
with charged grains. It only has a significant effect for pos- 
itively charged grains, for which the coulomb repulsion im- 
plies an extremely small rate of collisions with ions anyway. 
We still account for the electric dipole moment for the sake 
of completeness. 

5.5 Collisions with ions, neutral grain 

In that case the Coulomb potential vanishes, and the "image 
charge " potential has to be taken into account. We carry the 
calculation using the same assumptions as in the previous 
section : slowly rotating spherical grain, with radius acx- 
Taking /x as the polar axis for spherical polar coordinates, 
the interaction potential of the ion in the dipole and induced 
dipole field of the grain is given by 



v{r,e) = 



y2 2 3 



+ 



ZiQcji cos 6 



(94) 



2r2(r2-aL) 

The considerations that lead to the third constant of motion 
A hold again. The energy of the trajectory is 



where Vc«{r) is given by 



(95) 



Kff(r) 



2r^{r 



+ 



■m{bvY + 2ZiqclJ.cost 
2^ 



(96) 



Following DL98b, we define the dimensionless parameter 
^2 _ 2Zfq! 



Clcx 

kT' 



(97) 
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which describes whether the image charge attraction domi- 
nates over the thermal energy ((f>^ 1) or the thermal energy 
dominates {4> <Si 1). The effective potential can be written 



Kff(r) = kT 



J,2 4 



+ (U C + ^ cos 6*00) — 



(98) 



where fl , c, and u were defined in equations (|82|) and (|83|) . 
A study of this potential leads to the following condition for 
collision : 



max 

(c, u) 



~-l 



(99) 



The collision and excitation rates are obtained as in Eq. 
H86p . (|87p and (|88|) . One can then obtain the normalized 
damping and excitation rates for collisions of ions with a 
neutral grain: 



(100) 
(101) 
(102) 





= 0) 




= 0) 




= 0) 




I 4fl 



erf uo 



4 12 4 



(103) 



2/i0^ + (/)(2/i - 7)uo 
16^^ ' 



^ - 12/i + 15 + 20^) erf mo (104) 



(105) 



Note that in the limit /i ^ we recover DL98b result, as 

2 

-2 



0) 
0) 



2 ^ 

Q^l/2 



' + Y + c?(m 



(106) 
(107) 



However, the parameter fl is not small in general, as we 
saw in Eq. (|93p . so the effect of the dipole moment on the 
trajectory cannot be neglected. Note that we also have 
« 18 aZ^^^T^^^^ ■ The net effect of the dipole moment is 
to increase the collision and excitation rates, as can be seen 
from Fig. |4] In contrast to the case of charged grains, the 
electric dipole moment does have a significant effect and 
cannot be discarded. 

The effect of the dipole moment is always to increase 
both the collision and excitation rates, for both charged and 
neutral grains. This can be understood as follows. When the 
dipole moment vanishes, ions with a given velocity at infinity 
V collide with the grain if their impact parameter is such that 
b < 6max(i'). The effect of the dipole moment is to make a 
smooth transition from non-colliding to colliding trajectories 
: all ions with impact parameter b < bi{v) collide with the 
grain, a fraction {X{b,v) + l)/2 of those for which bi{v) < 



b < b2{v) do collide, and none of the ions with b > b2{v) 
collide, bi and 62 are such that 61 < &max < &2- As a result, a 
fraction of trajectories for which bi{v) < b < bmaxiv) do not 
lead to collision anymore (compared to the vanishing dipole 
case), and a fraction of trajectories 6max(w) <b < b2{v) now 
lead to collision. The suppressed colliding trajectories have 
a lower rate of collision and angular momentum than the 
added colliding trajectories. Thus the net effect of the dipole 
moment is to increase the collision and rotational excitation 
rates. 



6 PLASMA DRAG 

DL98b computed the effect of torques from passing ions 
on the electric dipole moments of the dust grains, which 
they named "plasma drag". They computed this effect for 
straight-line trajectories (the "Born approximation"). Here 
we include the full hyperbolic trajectory in the case of 
charged grains. We also account for the rotation of the grain 
explicitly. Moreover, we do not include trajectories leading 
to collisions, as they will give away their entire angular mo- 
mentum through collision, which we already accounted for. 
A precise calculation is important because plasma drag is 
one of the major excitation processes in some environments. 
Note that treatments of the plasma dra g effect that treat 
the plasma as a linear dielectric medium ('Ragot"2002') with 
the drag due to the imaginary part of the dielectric constant 
'it{u}, k) implicitly assume the Born approximation and do 
not capture the effects considered here. 

We will find that the straight-line approximation usu- 
ally overestimates the plasma drag. In the case of posi- 
tively charged grains, there is a range of impact parameters 
where the ion trajectory is deflected away from the grain, 
thereby suppressing angular momentum transfer. For nega- 
tively ch arged grains, ions can be foc used by electrostatic at- 
traction. I Anderson fc WatsonI l| 19931 ) argued that this is not 
a significant correction because the increased torque during 
close approach balances the shorter interaction time since 
the ion gains kinetic energy as it is attracted to the grain; 
however, we will see that in these cases there is a cancella- 
tion of angular momentum transfer in different parts of the 
trajectory that leads to reduced drag. For very special cases, 
the grain can corotate with the ion during close approach 
leading to an enhancement of the plasma drag, but this oc- 
curs for only a narrow range of impact parameters and does 
not compensate for the reduction of plasma drag that we 
find in other regimes. 

As in DL98b, we find it easiest to directly compute 
the plasma excitation Gp and use the fluctuation-dissipation 
theorem to infer R, = Gj,. 



6.1 Charged grain 

We consider the trajectories of positively charged ion (charge 
Zi > 0) in the electric potential of a charged dust grain 
(charge Zg 7^ 0). The trajectories are not strictly hyper- 
bolic due to the presence of the electric dipole potential (see 
Section I5.4|l . However, we saw that it has little influence 
on collisions and we will neglect its effect on the trajectory 
here, assuming they are hyperbolic and determined by the 
Coulomb potential only. The eccentricity of the hyperbolic 
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Figure 4. Fi{Zg = 0) (left panel) and Gi{Zg = 0) (right panel) for several values of the electric dipole moment, in CNM conditions, 
Eq. llT62l l. 



trajectory of the ion will be denoted e (as opposed to the 
elementary charge ^e)- 

Let the ion trajectory (a hyperbola) be in the {e.y,ez) 
plane, symmetric about the ey axis. The ion position is given 
by 



r = rer = r(a) (cos aey + sin a ez) . 



(108) 



The hyperbolic trajectory of impact parameter h and veloc- 
ity at infinity i; can be described in polar coordinates as 



r{a) 



e cos a 
P 



a e {—cte, Otc 
ce G (Qe, 2tv — c 



(Zg > 0), 

(Zg < 0); (109) 



the eccentricity and semilatus rectuirlfl of the trajectory are 



and 



= b\/e 



(110) 



The range of longitudes a of the trajectory are determined 
by the limiting case 



1 

arccos — . 

e 



(111) 



The longitude can be related to the true anomaly / familiar 
from planetary dynamics by q = / for repulsive {Zg > 0) 
cases and a — tt + f for attractive {Zg < 0) cases. We 
will need the following expression for the time t{a), valid in 
both cases (for the case of a n attractive potential, see e.g. 
iGevhng fc WestermanI (Il97ll ). Eq. (2.4.12)): 



t{a) 



1 



V e + 1 



e + 1 

e - 1 



In 



tan ■ 



tan ■ 



2e tan ■ 



tan^ f 



e-l 
e + 1 



(112) 



* The positive and negative cases of Eq. I I109I I could have been 
unified by taking the negative branches of the square root in 
Eq. miOI I. however we have not taken this route. 



In order to characterize the torque on the grain, we must 
first take the unit vector in the direction of grain rotation. 



sin 6 sin Sky + cos ( 



(113) 

so that — uiBui. We use {9, cp) to parameterize the (general) 
direction of rotation. We define the other two axes: 



cos y cos < 



cos y sm (p By 



and 



— — sm(l> Ex + cos(j) By. (114) 

In this system the electric dipole moment of the grain is 

fj, = + fi± [cos{ujt + x)e-g + sin{ujt + x)e<f>] , (115) 

where t = is taken to be the time when the ion is at the 
closest approach (i.e. r || By) and x ^ [0, 27r) is the random 
angle that /ix makes with Be at that time. 

The ion electric field exerts a torque on the grain dipole 
moment: 



7— = ax E 
dt ^ 



.ZiQc 



-/I X Br 



(116) 



-fJ, X Br- 



(117) 



Using the conservation of angular momentum, r a = bv, we 
can rewrite: 
div _ Ziq 
da Ibv 

We project that along the direction of b^ 
da 



■ ^ cos(aji + x) cos a cos (f> 



Ibv 

' &m.(u)t + x) (cos a cos Q sin < 



sm a sm 



0)](118) 



Expanding the sines and cosines, we integrate over the tra- 
jectory. We keep only the parts of the integral for which the 
inbound and outbound parts do not cancel, i.e. those which 
are even under a — > —a (Zg > 0) or a ^ 27r — a (Zg < 0); 
note that t{a) is even. We are then left with 



5uj« = 



(sin X cos 9 sin — cos x cos < 



Ibv 

X / cos (jji; COS ada 



cos X sin ^ / sinujtsmada 



(119) 
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In order to find the plasma excitation coefficient, we 
need to sum over collisions. We begin by averaging Jojy 
over solid angles for {9, cji) and over angles for x- The result 
is 

2 

(120) 



3 V Ibv 
We have defined the integral 



I I — ,e,Zg 

V 



cos Lut COS ada j 



+ sin cut sin ada^ , 



(121) 



where the integration limits are given by < a < cte (Zg > 
0) or Oe < a < TT (Zg < 0). Note that I only integrates 
over the inbound part of the trajectory; the outbound part 
is equal by symmetry. 

The excitation rate due to plasma drag is then given 

by: 
dAojif 



dt 



oo roo 



dv / 27r& db rii Anv 

■'i)max(l') 



Ihv 



. uib 

xX I e, Zg j , 



(122) 



where 6max(i'), the maximum impact parameter for collision 
to occur, is defined as 

OmaxW-| a,,^l^{2kT/mv2)i, mv'' /2kT > ^ ^ 
where ^ — ZgZiq^ / a^xkT . 

Note that technically the integration over impact parameters 
should stop at the Debye length 



kT 



: 398 



n<p/0.03cm-3 



1/2 



(124) 



However, we will see below that the integrand vanishes ex- 
ponentially for 



b> V ui 4.5 x 10 a_„ , / cm 

V rrii Vth oj 



(125) 



which is much smaller than the Debye length. 
Converting this into an excitation coefficient, we find 



Gv^ — ^l I . I X gp I \l uj I (126) 



where 

g,{^,n) = 2^e-" du j^^^^^ ^l(^l^,e,Zg ) , (127) 
where the eccentricity is given by 



(128) 



Note that we recover DL98b resullQ in the limit J = 1. 

^ DL98b include a term due to the parallel component of /x which 
is not relevant as it only leads to excitation perpendicular to lj. 



This expression has to be averaged over the grain charge 
and summed over all present ions. 

Straight line limit for I 

In the limit e oo, it is easier to express the integrals as a 
function of time, using 

y b 



V t 



sm a = 



da — 



1 



and 



dt. 



1 + 



(129) 



In this case, the first integral for I reduces to 



cos uit cos a da — — 



cos ujt dt 
[1 + (i;t/6)2]3/2 



uib ( cub 

V \ V 



(130) 



where Ki is a modified Bess el function of the second kind. 
[Here we used Eq. (9.6.25) of lAbramowitz fc Steguiil (Il972l ) 
with 1/ = 1, z = 1, and x — u;b/v.] The other integral is 



/ 



sin ujt sin a da 



vt 



sin Lot 



b [1 + (i;t/6)2]3/2 
r sin XT dr 



dt 



,0 (1 + ^2)3/2' (^31) 

where x = u)b/v. Since r(l + t-2)-3/2 jg ^jjg derivative of 
— (1 + r^)"^''^, we can integrate by parts and find 



+ 

Jo 



X COS XT dr 



(132) 



The boundary terms evaluate to zero, and the second inte- 
gral can again be evaluated t o xKo{x) using Eq. (9.6.25) of 
lAbramowitz fc Steguiil l| 19721 ) with u = and z ^ 1. Thus 
we have 



x^lK^oi^) + K^ix)], 



Lob 

V 



Note that when ui 
' u)b 



we recover DL98b result, i.e. 



= 0, e 



Zg ] = 1. 



(133) 



(134) 



We moreover have an exact functional shape for the cutoff 
at large rotation rates. 

Non rotating grain limit for T 
It is straightforward to show that 
1 



T ( — = 0,e 

V 



(135) 



for both positively and negatively charged grains. Thus, the 
nearly parabolic trajectories e — 1 <C 1 are suppressed by a 
factor ~ 2(e - 1). 

The numerical calculation of X in the general case is 
tricky because it involves integrating an oscillating function 
which frequency goes to infinity at one limit of the integral, 
as t(a Qe) ^ oo. We refer the reader to Appendix |3 
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for the description of the method used for numerical com- 
putation. Fig. [5] shows the resulting dimensionless torques. 
An important feature is that for negatively charged grains, 
ions with nearly parabolic trajectories may corotate with 
the grain which results in an enhanced torque. 



6.2 Neutral grain 

The exact calculation of the trajectory in the electric dipole 
potential and the "image charge potential" is untractable 
analytically, and would require a heavy numerical calcula- 
tion. Therefore, we will make the following simplifications. 
First, we neglect the effect of the electric dipole moment on 
the trajectory. This assumption is somewhat cavalier, as we 
saw previously that the electric dipole moment may signifi- 
cantly afi'ect the ion trajectory in the case of a neutral grain. 
Furthermore, although trajectories in the "image charge po- 
tential" will be curved in general, we will consider them to 
be straight lines. Thus, we will approximate the torque given 
to the grain by Eq. (fMJ)) . where I is given by Eq. (fT^ . Col- 
liding trajectories should not be taken into account for the 
plasma drag excitation rate. Thus, we integrate the torque 
only over trajectories with impact parameter b > bmax(t'), 
with 



1+- 

U 



(136) 



(see DL98b Eq. (B24) and the definition of <^ Eq. (|97)l 'l 

Therefore, in the case of neutral grains, we have 
rii I rrii f ZiqcfJ._ 



Gp{Z, = 0) 

where 

9p (0, ^) = 



2kT 



uj (137) 



2ue " Au 



(138) 



The normalized excitation rate for plasma drag for Cold 
Neutral Medium conditions, Eq. (|162|l , is shown in Fig. |6] 



7 INFRARED EMISSION 

A dust grain absorbs light and reemits it in the infrared. A 
rotating grain will also radiate angular momentum, which 
damps its rotation. 

DL98b compute this damping rate by modelling the 
grain as composed of six rotating dipoles. We give here a 
more accurate calculation, using the correlation functions of 
the dipole moment in the grain frame. Our result is a factor 
of two greater than that of DL98b. We present a classical 
calculation in this section; the quantum calculation is pre- 
sented in Appendix |B] and gives the same result. 

The rates of emission of energy and of angular momen- 
tum by a varying electric dipole moment: 

2 



E 



2 ..2 

3^^ 



and L 



3c 



-p X p. 



(139) 



We denote the coordinates of the dipole moment in the frame 
corotating with the grain with unprimed indices, and the 
ones in the "lab frame" with primed indices. Take a grain 
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Figure 6. Normalized excitation rate due to plasma drag Gp(a;th) 
for a neutral grain, a positively charged grain, and a negatively 
charged grain in CNM conditions (Eq. I|162|l ). evaluated at the 
"thermal rotation rate" ojth = {2kT / 1)^^'^ . All of them are lower 
than estimated by DL98b. It is clear that the positively charged 
grains are much less excited than the neutral and negatively 
charged grains. The kink at 6A is due to the change of grain 
shape. 



rotating around the z-axis, without precession, with angular 
frequency u. We have 

Px' ~ cos ut px — sin ut Py , 

Py — sin uit Px + cos Ljt Py, and 



(140) 



A straightforward calculation leads to the following expres- 
sions in the lab frame : 

= pl + pI + pi + 'iuj [P^Py - PyP^) 
W \4.{pl+pl) 



2 {pxPx +PyPy)\ 
+AUJ^ {PxPy - PyPx) + Uj'^ {pl + pI) 



and 

(p X p)^ 



PxPy 



■ PyPx +U}[2 (pl + pI) 



+ ?,U>'^ {PxPy - PyPx) + Uj'"^ {pl + pI 



(141) 

(PxPx +PyPy)] 
(142) 



Since we are interested in the statistical properties of 
the emission, we define the unequal-time dipole moment cor- 
relation function in grain coordinates. 



a,{r) = {{p,{t) - {p.))(pAt + t)- {p,})) , 



(143) 



where (pi) = is just the constant dipole moment of the 
grain. We further assume statistical spherical symmetry of 
the dipole moment in the grain coordinates, i.e. dj = C'Sij. 
(For a planar grain, the values of the correlation functions 
depend on the in-plane or out-of-plane character of the vi- 
brational modes and may be anisotropic. However if the in- 
frared emission arises during thermal spikes when the grain 
is not rotating around its axis of greatest angular momen- 
tum, we expect the isotropic analysis to be a good approx- 
imation.) The average values of the previous formulae then 
becomqlj 

Expectation values of derivatives such as (p^) can be expressed 
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Figure 5. Contour levels of e, Zg > 0) (left panel) and , e, Zg < 0) (right panel). Both show that I goes to unity for slowly 

rotating grain, straightline trajectories, and vanishes for rapidly rotating grains or nearly parabolic trajectory. In the case of negatively 
charged grains, though, there is a visible corotation regime, where e — 1 <g 1 and uib/v{e — 1) ~ 1, for which the ion and the grain 
approximately corotate, enhancing the torque given to the grain. 



(p ) = 3C""(0) - 12lu^C"{0) + 2lu''C{0) (144) 
and 

{P X p), = ~6ujC"{0) + 2uj'-^C{0), (145) 

where ' denotes the derivative of the correlation function 
with respect to r. 

The Wiener-Khintchine Theorem relates the corre- 
lation functions to the Spectral density S,^, C{t) = 

Sv cos(2-K V t) Au . Plugging back into Eqs. (1144 |l and 
(|145p . we get 



roc 

{p'^) = [3{2Tnyf + 12lo'^ (2niyf + 2lo^] du (146) 
Jo 

and 

{p y-p)z= I [&u}{2TTvf + 2lu^] Av. (147) 
Jo 

Now making use of the assumption that the grain rotates 
slowly, i.e. that Vrot = lu/2ti <^ = typical frequency of 
emission, in the infrared, we get, at the lowest order in 
I'rot/vQ, the average total power and average total rate of 
radiation of angular momentum; 



d£;\ 2 2 



At 
and 
'AL. 
At 



3c3 



3c3 



{2-KufS, Av 



{P X p). 



4uj 



{2-KvfS^ Av. 



(148) 



(149) 



If one knows the infrared power radiated per steradian 
per frequency interval Fu, such that 



AE 

At 



■ 4tt 



F^Au, 



(150) 



one can deduce the rate of angular momentum loss through 
infrared emission: 



in terms of correlation functions via integration by parts. In this 
case, (p^) = {(pxpx)) — (Pxpx)- The first term vanishes for a 
stationary process, and the second is — C"(0). 



AL, 
At 



2lo 



F^ 



Af. 



(151) 



This result is twice as big as the one given in DL98b. [The 
difference occurs because DL98b modeled the dipole fluctu- 
ations with six uncorrelated rotating dipoles, one rotating 
each direction in the xy, yz, and xz planes. They counted 
the radiated power from all six of these, but only considered 
the angular momentum loss from two of them (in the xy 
plane). The dipoles rotating in the xz and yz planes con- 
taining the rotation axis also emit net angular momentum 
however, and if they are considered one recovers the factor 
of 2.] 

This classical treatment does not predict the rotational 
excitation from the recoil given by individual photons, which 
is a quantum effect. As in DL98b, we set 



dAL" 
At 



AN, 



phot 



At 



F. 



Aiy. 



(152) 



The normalized damping and excitation rates are then 



-Fir 

GiR, 



2th 

TT-f Jo 

h m 
6ttI kT jQ 



Fv 

— ^ dv and 



F, 



Av. 



(153) 



We calculate the infrared emissivity of PAHs and small car- 
bonaceous grains using the "thermal continuous" approx- 
imation, described in DLOl. They indeed show that this 
treatment leads to spectra very close to those predicted by 
the exact statistical treatment, and has the advantage of 
being computationally much faster. We obtain the steady- 
state energy distribution function and then get the infrared 
emissivity, as explained in DLOl. 

We checked numerically that we recover DL98 result for 
low values of the radiation field intensity : Jir, Gir oc %. 
However, their result for high values of the radiation field 
(Fir oc x^''^! Gir oc X"^^) relies on the fact that the ab- 
sorption efficiency Qi, oc v'^ at the characteristic frequencies 
of infrared emission. This is not valid anymore for high ra- 
diation fields, which offset the emission spectrum to higher 
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frequencies, where the absorption efficiency has not a simple 
dependence on frequency anymore. We show the resulting 
infrared emission and damping coefficients in Fig. [T] 



8 PHOTOELECTRIC EMISSION 

An electron ejected from the grain carries away an angular 
momentum along the rotation axis (z-axis) equal to : 



ALz = ruep (w^ - puj) 



(154) 



where v'^ is its tangential velocity in the grain frame. From 
this we deduce that 



J nt 



mn 2Tva^nH^/2kT/TvmH ' 



(155) 



where Jpe is the photoemission rate and was described in 
section 13.41 The excitation rate can be obtained by first 
noticing that the rotational velocity is much smaller than 
the velocity of ejected electrons: 



au) <ti 

so that we have, up to small corrections 



A r2 2 2 /2 



(156) 



(157) 



We assume a cosine-law directional distribution for escaping 
electrons, so that {v'^) — \v1, where we denote Ve the aver- 
age velocity of the electron at the grain surface. The latter 
satisfies 



{Z, + l)ql 



En 



(158) 



where Epe is the average energy at infinity of the photoe- 
jected electron. We finally get 



/ A r 2\ 2 ^ "ca; J- 2 '"-e ilea; 



So the normalized excitation rate is 



Gpe — 



4nH (STTmnkT)^^^ alkT 



{z, + i)d 



(159) 



, (160) 



where Fpe is the heating rate due to photoemission of elec- 
trons, obtained from WDOlb. 



9 RANDOM H2 FORMATION 

DL98b showed that the random formation of H2 molecules 
on the grain surface does not make a major contribution to 
rotational excitation. We use their result: 
(J(J + l))ft2- 



(161) 



2mYiEfa% 

where 7 is the efficiency of H2 formation, y = Irin^/nn, 
Ef ~ 0.2eV is the average translational kinetic energy of the 
nascent H2, and (J(J -I- 1)) ~ 10^ gives its average angular 
momentum. 



10 RESULTING EMISSIVITY AND EFFECT 
OF VARIOUS PARAMETERS 

Throughout this section and unless otherwise stated, we will 
take as a fiducial environment the CNM parameters specified 
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Figure 8. Rotational distribution function for a grain radius a = 
7 A, in CNM conditions, for a single value of the dipole moment 
/ii(10-''cm) = 9.3 debye. The plot compares the solution of the 
Fokker-Planck equation with the DL98b Maxwellian approxima- 
tion (DL98b Eq. (57) used with our F, G). Note that DL98b pre- 
scription gives (< ijj^ >)^^^ = 2ttX 10.7 GHz, which is in excellent 
agreement with the value we get, (< i^;^ >)^^^ = 2it X 10.9 GHz. 
However, the shape of the distribution function is significantly 
diff'erent. 



by 
riH 

X 



30 cm"'' , T = 100 K 
n(H+)/nH = 10"=^ , xc 
u/uiSY^¥ = 1 , 7 = 0. 



; n(C+)/nH = 3 X 10"* 
(162) 



We also take the rms intrinsic dipole moment to be 
{^J.'i)^^^ [a — 10"^cm) = 9.3 debye. 



(163) 



For the size distribution parameters, we use those given by 
WDOla for a ratio of visual extinction to reddening Rv = 
3.1, and a carbon abundance in the log-normal distributions 
&c = 6 X 10"^ 

This section is intended to give some intuition into the 
effect of various parameters on the spinning dust spectrum. 
However, the reader should keep in mind that environment 
parameter space is many-dimensional, and changing several 
parameters at once may lead to modifications that are not 
superpositions of the effects described here. 



10.1 General shape of the rotational distribution 
function 

The rotational distribution function is obtained as described 
in section |4] We remind the reader that the Fokker-Planck 
equation is not stricltly valid for the smallest grains, for 
which impulsive torques are important. It however still de- 
scribes their rotational distribution function with more accu- 
racy than a simple Maxwellian. Moreover, DL98 showed that 
impulsive torques may be neglected for grain radii a ^ 7A. 
In Fig. |8] we show that the rotational distribution function 
obtained by the Fokker-Planck equation differs significantly 
from a Maxwellian. It has a sharper cutoff at high frequen- 
cies due to the proper accounting for rotational damping 
through electric dipole radiation. 
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a (cm) a (cm) 

Figure 7. Infrared omission damping and excitation coefficients -Fir, Gir, in CNM conditions (Eq. I|162|l '). compared with DL98b result. 
The difference is mainly due to differences in grain absorption efficiencies and the calculation of the infrared spectrum. We used Li & 
Draine (2001) absorption efficiencies and DLOl model to compute the infrared emissivity. The kink at 50 A in our result is due to a 
change in optical properties of dust grains. The kink around 50 A in DL98b result is due to the change from constant temperature limit 
(larger grains) to thermal spikes limit (smaller grains). The fact that they coincide is purely chance, and would not be necessarily the 
case for other environment conditions. The discontinuity at 6A is due to the change in grain shape. 



In what follows we will analyze the effect of various pa- 
rameters on the spinning dust emissivity. As can be seen 
from equation (|33p and the expressions derived next for the 
normalized damping and excitation rates, the rotational dis- 
tribution function has complex dependencies on all grain 
and environment parameters. To get some intuition on the 
physics of spinning dust and the influence of each parame- 
ter, we will rely on a simplified expression for the rotational 
distribution function in the following sections: 



/„(..) ocexpl-- — 



TH 



3Gred \2kT) 



(164) 



where we approximate the plasma drag excitation rate 
(which is in principle a function of uj) by the constant 



1/2 



(165) 



In our analysis we will also neglect the charge-displacement 
induced dipole moment as it has a minor contribution. Of 
course the actual rotational distribution function and emis- 
sivity are computed using the exact equations developed in 
this paper. 

For a given grain radius a and intrinsic electric dipole 
moment pi, the power radiated is (a; /Xi) oc fa{2-KV\ /li). 
It is straightforward, from Eq. (|164|l , to show that the peak 
frequency is given by 



t^pcak 



G 



1/2 



6kT 



where we defined the parameter 



^ ~ l?2 



8G TH 

Ted 



(166) 



(167) 



which denotes the non-Maxwellian character of the distri- 
bution function. 

For ^ ^ 1, the distribution is nearly Maxwellian, 



/aHocexp(^-|^), 

and the peak frequency is given by 



gV^' 1 



6kT 



(168) 



(169) 



Moreover, the total power emitted by a single grain 
ja (X fj,^ J Lo^ fa{io)duj has the following dependence: 



i,(XM'(fj (170) 
For ^ 3> 1 the distribution is strongly non-Maxwelliarf"1. 



fa{L0)(X exp 



3Gred \2kTJ 



and the peak frequency is given by 



i^pcak ' 



V 2-rH / 2-K 



(C » 1). 



The total power is then given by 

2 Grod li- ^ i\ 

JaOCfl T (C>1)- 



(171) 



(172) 



(173) 



In Fig. |9]we show the rms rotation rate < u'^ >i/2 
a function of grain radius. As can be expected, the small- 
est grains are rotating with the greatest angular velocity, 
as they have the smallest moment of inertia. Consequently, 
they radiate at the highest frequencies, and constitute the 
peak of the spectrum. Therefore, we will use Eqs. (|164p to 
(|173p for a grain of radius ttmin = 3.5 A to evaluate the effect 
of various parameters on the emissivity. 

We finally remind the reader with the dependencies of 
characteristic timescales : 



Interestingly, lEricksoiil l|l953) had already obtained a result 
similar to Eq. II171II with a Fokker-Planck equation. 
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Figure 9. rms rotation rate uiims =< ^^'^ >i/2 as a function 
of grain radius a, for CNM conditions. The rotation rate the 
grain would achieve if it were rotating thermally (in that case 
i.^rms = \/3kT/I) is also shown. It can be seen that the grains 
rotate subthermally. The kink at 6 A is due to the change in grain 
geometry. 



Figure 10. Power radiated by one grain of radius a = 3.5 A 
in CNM conditions. The dotted and dashed lines show the 
contributions of various values of the intrinsic dipole moment, 
which is assumed to have a gaussian distribution with rms value 
fj,0 =< fJ-l >^^^ (a). The solid line is the total power. 



TH oc nn'T-'/^ , Ted oc /i-'r-' (174) 

10.2 Emissivity 

Once the rotational distribution function is known, as a 
function of the intrinsic electric dipole moment, fa{uj;jJ,i), 
one can get the power radiated by a grain of radius a by 
averaging over the intrinsic dipole moments gaussian distri- 
bution P{fii) defined in equation ([8]). One gets : 

/2 2 6 
d^i^Pip,) - 2tt fa{io;fi,) (175) 

where IJ,\ = f for spherical grains, and fi^ = fi^ for cylin- 
drical grains. 

The overall effect of averaging over the dipole moments 
distribution is to broaden the spectrum, as can be seen in 
Fig. IIUI The peak frequency remains approximately equal 
to that of Pi,{ni = (fMi)^^^). We will discuss the effect of the 
rms intrinsic dipole moment in section flU. 31 

The emissivity per H atom is then obtained by inte- 
grating the power radiated by each grain over the grain size 
distribution function, described in section The emissiv- 
ity for the CNM environment is shown in figure 1111 Note 
that the grain size distribution directly weights the spec- 
trum, and thus needs to be known with accuracy, which is 
not quite the case yet for the very small grains. 

10.3 Effect of the rms intrinsic dipole moment 

Varying the rms intrinsic dipole moment affects the spec- 
trum in three main ways. First, it affects the total power 
radiated, as oc fi^. Then, it affects the non-Maxwellian 
character of the distribution function, as Ted oc fi^. Finally, 
it affects the rotational damping and excitation rates essen- 
tially through plasma drag, which has Gp oc fi^ (the effect 
on Gi is not as important). 
Throughout the range of values considered. 




V (GHz) 

Figure 11. Spinning dust emissivity for CNM environment. 
Contributions from grains of various sizes are shown (a_7 = 
a/(10~^ cm)). The grain size distribution parameters are taken 
from WDOla with Ry = 3.1 and 6c = 6 X 10^^. 

1 debye < /i^lO"^ cm) < 100 debye, (176) 

and in CNM conditions, the distribution function remains 
strongly non-Maxwellian : ^ > 60. Therefore, we can use 
the strongly non-Maxwellian limit Eq. (I172|l to evaluate the 
peak frequency. 

Low dipole moment limit 

For low values of the electric dipole moment, plasma drag 
has little effect on both the rotational damping and excita- 
tion. Therefore, F and G are roughly independent of fi, and, 
from Eqs. (I172|l . (|173|l and Ted oc we get 

I'poak (m» ^ 0) oc ^7^/'' (177) 

j/nu ilJ-i — > 0) constant (178) 

One can see in Fig.[T2]that Eq. (|177|l is quite accurately sat- 
isfied. The total power has a weak dependence on /Xi for low 
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values of the intrinsic dipole moment, but is not strictly inde- 
pendent of it, which comes from the multiple approximations 
made in this analysis (neglecting the charge displacement- 
induced dipole moment, and using Eq. (|173|) for the total 
power, after integration over the size distribution, instead of 
the total power radiated by a single grain). 

High dipole moment limit 

For high values of the electric dipole moment, plasma drag 
dominates both rotational damping and excitation. There- 
fore, G ~ Gp oc fi^ so we get 



{fii oo) ^ constant 
j/nn (pi, oo) oc , 



(179) 
(180) 



which describe approximately the behavior observed in 



10.4 Effect of the number density nn 

The main effects of the number density are : 

• Changing the relative contribution of gas-induced 
and radiation-induced rotational damping and excitation. 
For very low number densities. Fir and Gir oc th oc rijj^ 
dominate over other rotational damping and excitation 
rates. For high densities, plasma drag and collisions are 
dominant. Note that the charge distribution is also modified 
as the higher the density, the more important is coUisional 
charging compared to photoemission. As a consequence, the 
grains are positively charged at low densities, and tend to 
be negatively charged at high densities, due to the higher 
rate of collisions with electrons. 

• Influencing the non-Majcwellian character of the rota- 
tional distribution function. The higher the number density, 
the closer is the distribution function to a Maxwellian. Nu- 
merical calculation shows that starting from CNM condi- 
tions and varying only nu, we transition to the Maxwellian 
regime (C < 1) if nn > 10^ cm'^ 

Low density limit 
For very low number densities, the distribution is highly non- 
Maxwellian and we can use Eqs. (|172|l and (|173p . with G — 
Gir, to estimate the peak frequency and total power. As 
Gir/th is independent of nn, both the number density and 
total power should asymptote to a constant value. We can 
estimate numerically the peak frequency in CNM conditions 
and get: 

i^peak {nn ^ 0) ^ 13 GHz (181) 
j/nu (nn ^ 0) ^ constant (182) 
which is in good agreement with Fig. 1131 

Intermediate densities 

Over the range 10^ cm~^ < tih < 10^ cm~^, gas pro- 
cesses are dominant over infrared emission, so F, G are 
roughly independent of tih. In addition, the distribution is 
still strongly non-Maxwellian. Using Eqs. (|172|) and (|173|) . 
we thus find 



The kink around tih ~ 3 x 10'^ cm'^ is due to our discontinous 
treatment of the evaporation temperature for high densities 
(see section [5.1.4p . and to the replacement of the integration 
over all grain radii by a discrete summation when numeri- 
cally computing the spectrum. Therefore the spectra should 
not be considered as very accurate in that region. 

High density limit 
For very high number densities, the excitation and damp- 
ing is dominated by gas processes, and the electric dipole 
damping becomes negligibly small, so that the rotational 
distribution function is actually a Maxwellian, although not 
thermal. Using Eqs. (fTMIl and (IT70|l . we find 

I'peak (nn ^ oo) ^ 150 GHz (185) 

j/nu ("-H oo) ^ constant (186) 

10.5 Effect of the gas temperature T 

Temperature has a less obvious effect on the spectrum and 
we need to analyse in detail every damping and excitation 
process. It turns out the charge distribution of the smallest 
grains varies very little over the range of temperature con- 
sidered 1 K < T < lO'^ K and they remain mostly neutral 
throughout this interval. The distribution remains strongly 
non-Maxwellian for T greater than a few K. 

Low temperature limit 

At very low temperatures, the dominant excitation process is 
collisions with ions. Indeed, the grains being mostly neutral, 
the ions interact strongly with the electric dipole potential. 
As /i oc and (p oc T~^^^, one can see from Ea. (|102|l and 
(|104p that Gi oc T~^. Plasma drag has also Gp oc in 
principle but this becomes a shallower power law at low tem- 
peratures as the interaction timescale becomes longer than 
the rotation timescale. We find numerically, though, that 
roue hly G oc as G is not strictly equal to Gi (colli- 

sions with neutrals are also significant at low temperatures) . 
Using Eqs. (fT69)) . (tTTOj) . and (IT74)| we find 

i^poak (3 K < T < 10^ K) 35 GHz (187) 

j/nn (3 K < T < 10^ K) « constant (188) 

Note that for extremely low temperatures, the distribution 
would become Maxwellian, and one would get, according to 
Eqs. ((T721), <fT75)) . 



I'poak (T 0) oc r'/' 

j/nu (T -» 0) oc 



(189) 
(190) 



3^ 1/4 
oc n-J 



i^pcak (10 cm < riH < 10 cm ) 
j/nn (10^ cm~"^ < nn < lO** cm~"^) oc nn 



(183) 
(184) 



which can be guessed at the extreme low temperature end 
of Fig. [131 Temperatures below ~ 3 K are of course unphys- 
ical, but for other environmental conditions than those of 
Eq. (|162p . the behaviour discussed above could take place 
at higher, observed temperatures. 

High temperature limit 

At very high temperatures, collisions with neutrals are the 
dominant damping and excitation process. The CNM envi- 
ronment being mostly neutral, _F„ — ^ 1 and G„ —> 1/2 at 
high temperatures (Gi"^"' oc T^v/T 0). Moreover, the dis- 
tribution becomes strongly non-Maxwellian, as ^ oc T^^^. 
We therefore obtain 



i'poak(r - 



200 T^^'^ GHz 



(191) 
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Figure 12. Effect of the intrinsic electric dipole moment on tlie peak frequency and the total spinning dust emission. Environment 
parameters are set to CNM conditions defined in equation (|162|l . Increasing the electric dipole moment decreases the peak frequency and 
increases the total power radiated. 



j/riH (T ^ oo) oc r^/^ (192) 

Fig.[l3]shows that these power laws describes the behavior ol 
the peak frequency and total power with very good accuracy. 

10.6 Effect of the radiation field intensity x 

The radiation field affects the spectrum through only two 
ways. First of all, it changes the charge distribution of the 
grains as an increased radiation implies a higher photoe- 
mission rate. Second of all, it affects the rate of damping 
and excitation through infrared emission, (and photoelec- 
tric emission, but this is subdominant) . 

Low radiation intensity limit 

In a low radiation field, -Fl-ad and Grad become negligible. 
The photoemision charging rate becomes insignificant com- 
pared with coUisional charging, and the charge distribu- 
tion function depends only on other environment parame- 
ters. Thus, one expects the spectrum to reach an asymp- 
totic shape for very low radiation fields. The distribution 
is strongly non-Maxwellian, and the dominant excitation 
mechanism is collisions with ions, whereas the dominant 
damping mechanisms are plasma drag and collisions with 
neutrals. Thus, we find 

i'peak(x^ 0) ^ 35 GHz (193) 

i/'^H (x ^ 0) ^ constant (194) 

The kink around x ~ 2 x 10~^ is due to our discontinuous 
treatment of the evaporation temperature for low intensities 
of the radiation field. 

Around x ~ 1 — 10, the grain becoming more and more 
positively charged, the collisions with ions start being less 
efficient, although still the dominant excitation mechanism. 
This results in a slight decrease in both fpeak and j/nu- 

High radiation intensity limit 

In a high radiation field, F ~ Fib., and G ~ Gir. Both Fm. 
and Gir are approximately (although not strictly) linear in 
X, as shown in DL98b for the thermal spikes limit (see their 
equations (31) and (44)). Thus, ^ ~ so the distribution 
becomes Maxwellian. The peak frequency and total emitted 



power are then given by Eqs. (|169p and (|170p . which imply 
that 

!^pcak(x ^ oo) « constant (195) 

j/i^H (x ~* oo) ~ constant (196) 

These asymptotic forms are not strictly valid because Fir 
and Gir are not strictly linear in x, and do not have a simple 
dependence on that parameter. 

10.7 Effect of the ionization fraction xu 

The Hydrogen ionization fraction affects the charge distri- 
bution by modifying the contribution from collisions with 
protons. It also changes the contribution of collisions with 
ions, neutrals and plasma drag. Characteristic timescales are 
left invariant, and ^ 3> 1 for any ionization fraction in oth- 
erwise CNM conditions. 

Low ionization fraction limit 

In that limit the rotational distribution function reaches an 
asymptotic form where collisions with protons and plasma 
drag due to protons can be neglected. However, there are 
still G^ ions in the gas so collisions with ions and plasma 
drag may still be important, although the dominant excita- 
tion process is collisions with neutrals. We find 

!^pcak (XH -^0)^ 30GHz (197) 

j/nn (xh — > 0) ^ constant (198) 

High ionization fraction limit 

In that case collisions with ions are the dominant excitation 

process. Using Eqs. (|172|) and (|173|l along with G ~ Gi oc 
sh, we find 

J^pcak(a;H ^ 1) ~ 90(s:h/0.1)'/''GHz (199) 

j/nn (xh ^ 1) oc sh (200) 

10.8 Concluding remarlcs 

We remind the reader that all the estimates in the previous 
section were given by assuming that the peak frequency of 
the spinning dust spectrum is determined by that of the 
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Figure 13. Effect of various environmental parameters on tlie peak frequency and the total spinning dust emission. When one parameter 
is varied, other environment parameters are set to CNM conditions defined in Eq. II162I I. See sections ll0.4l to[ l0.7l for a detailed description. 
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smallest grains, and that the total power follows the same 
dependence upon environmental parameters as the power 
emitted by the smallest grains. Therefore they should be 
taken as an aid to understand the physics of spinning dust, 
but not as an accurate description, which requires numerical 
calculations. 

The overall conclusion of this section is that varying a 
single environmental parameter may change the peak fre- 
quency by up to an order of magnitude, and the total emit- 
ted power by several orders of magnitude. There is therefore 
a very large range of possible peak frequencies and total pow- 
ers that can be produced by spinning dust radiation. Mul- 
tiphase environments, in particular, could emit very broad 
spinning dust spectra. Deducing the environment parame- 
ters from en observed spectrum could therefore be a difficult 
task. 

We show the spinning dust spectrum for various envi- 
ronments and compare them to DL98 results in Fig. 1141 



11 CONCLUSION 

We have presented a detailed analysis of the rotational exci- 
tation and damping of small carbonaceuous grains. We have 
refined DL98b results in the case of collisions, accounting 
properly for the centrifugal potential which increases the 
net damping rate. In the case of collisions with ions, we ac- 
counted for the effect of the electric dipole potential on the 
collision cross section. We found that this is a small effect 
in the case of charged grains, but that it may significantly 
increase the damping and excitation rates in the case of neu- 
tral grains. We evaluated the contribution of "plasma drag" 
by considering hyperbolic trajectories and rotating grains in 
the case of charged grains, and straightlines in the case of 
neutral grains. We corrected DL98b results for the damp- 
ing through infrared emission. Finally, we calculated the ro- 
tational distribution function by solving the Fokker-Planck 
equation. 

We believe our model provides a much more accurate 
description of the spinning dust spectrum than previous 
work. However, we would like to remind the reader of its 
uncertainties and limitations. First, our model only com- 
putes the total intensity of the emitted radiation and not 
the polarization, which would require an additional study 
of the alignment mechanisms for PAHs. Secondly, the dust 
grains properties are poorly known: 

• The size distribution and abundance of the smallest 
grains is uncertain, and in particular the nature of the cutoff 
at small grain sizes a ~ flmin can have a large effect on the 
spectrum. 

• The permanent electric dipole moments of dust grains 
are not directly constrained by other dust observables. Given 
that it cannot be computed from first principles, one may 
regard it as a free parameter (or parameters) of the spinning 
dust model. 

Thirdly, we made some simplified calculations in some cases, 
as an accurate calculation would have been intractable nu- 
merically or substantially complicated the code: 

• We used the Fokker-Planck approximation, which starts 
to break down for our smallest grains because a single colli- 
sion suffices to change the rotational state. We expect that 



the main consequence of a full treatment would be a tail 
in the emission spectrum extending to high frequencies, be- 
cause impulsive collisions would be able to increase the ro- 
tation velocities of the grains to > 2j^pcak before dissipative 
forces had time to act (an effect missed by the Fokker-Planck 
treatment). Therefore one should not place too much confi- 
dence in the many order-of-magnitude falloff at ~ 100 GHz 
seen in most of our models. (In many cases this will be unim- 
portant observationally since at high frequencies the vibra- 
tional dust contribution is dominant.) 

• In the plasma drag calculation, we neglected the electric 
dipole potential when evaluating the trajectory of ions, tak- 
ing the straight-line (neutral grain) or hyperbolic (charged 
grain) approximation. Relying on the study of collisions with 
ions, we may expect the dipole moment to have a small efi^ect 
in the case of a charged grain. On the other hand, its effect 
in the case of a neutral grain may be more important, as in 
that case the electric dipole potential provides the dominant 
interaction. 

• We assumed the evaporation temperature for the small- 
est grains was the "temperature" of the grain just after it has 
absorbed a UV photon. This is a physically motivated as- 
sumption but its validity is not established. The evaporation 
temperature can have a significant effect on the spectrum, 
as can be seen from Fig. [3] and one should be aware of the 
uncertainty in this parameter. Also, we assumed that colli- 
sions transition from being sticking to elastic, as the density 
exceeds a given threshold. Our model is therefore inaccurate 
in the transition region. 

• When calculating the infrared emission spectrum of the 
grains, we used DLOl "thermal continuous" approximation, 
which is not very accurate to describe the low energy part 
of the spectrum. Whereas these uncertainties are not impor- 
tant if one only wants the spectrum Fi, in the mid-infrared, 
they may lead to significant errors when calculating the cor- 
responding damping and excitation rates, which are propor- 
tional to / v~^Fvdv and / v~^Fvdv respectively. 

• We ignored systematic torques, although this may not 
be a major omission for the smallest dust grains. 

Despite these uncertainties, we believe that this model 
is the most complete thus far, and will be a useful tool for 
comparison to observations and testing the spinning dust 
hypothesis for anomalous microwave emission in various ISM 
phases. 
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Figure 14. Spinning dust spectra for several environment conditions : Cold Neutral Medium (CNM), Warm Neutral Medium (WNM), 
Warm Ionized Medium(WIM), Molecular Cloud (MC), Dark Cloud (DC) and Reflection Nebula (RN). The environments are defined in 
DL98b, Table 1. The thin solid line is the result of our calcu lation, the dotted line being DL98 prediction, and the dashed line is the 
free-free emission (the free-free gaunt factor were taken from ISutherlanj |l99i)). The parameters for the grain size distribution are : 
Rv = 3.1, fee = 6 X 10-5 for the diffuse CNM, WNM and WIM environments, and Ry = 5.5, be = 3x IQ-^ for the dense MC, DC and 
RN environments. The apparent systematic increase of power around the peak frequency for our result is mainly due to the grain size 
distribution we use, which has an increased number of small grains compared to that used by DL98 (compare e.g. Fig. 2 from WDOla 
with Fig. 8 from DL98b). Note that for the DC environment, for which rotational excitation and damping is dominated by collisions 
with neutral species (mainly H2 molecules), DL98 prediction largely overestimates the peak frequency and total power because they 
underestimate the damping rate (see Eq. II53I I and subsequent discussion). 
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APPENDIX A: PLASMA DRAG : NUMERICAL 
CALCULATION OF I IN THE GENERAL CASE 

The numerical calculation of T is tricky because it involves 
integrating an oscillating function which frequency goes to 
infinity at one limit of the integral, as t{a —> Oe) —> oo. 
Here we describe our implementation for both the positive 
and negative grain charges. 

Al Positively charged grains: Zg > 

We first make the change of variable 

= V7cot|, (Al) 

where 7 = (e — l)/(e + 1). The expression for the time is 
now: 



ujt{z) 



tub 



1 



V - 1 V -2-1 
The X-integral is then 
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T = 47 
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dz 



(A3) 



The functions inside the integrals are analytical on the 
complex plane, deprived from the branch cut [—1, 1] on the 
real axis and the two poles ±1^/7. The integrands are at 
least 0(2"^) as \z\ 00. Moreover, for y 0+, 



3? [iLjt(l - iy)] (X 3ff 
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-1 + ^V^ 

y J y. 
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(A4) 



Thus, using the fact that the integral over the lower right 
part of the complex plane vanish, we can replace our inte- 
grals by integrals over the axis 



< y < +00. 



(A5) 



Note that for e — » 1, T = 0{e — 1), as one may expect from 
almost parabolic trajectories if the grain repels the ion. Also, 
in the limit cub/v 0, T -> (e^ - l)/e'^. 

A2 Negatively charged grains: Zg < 

This time we make the change of variable 
1 a 

The expression for the time is now 
ujb 1 



(A6) 
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And we have: 
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(1 +722)2 

The functions inside the integrals are analytical on the com- 
plex plane, deprived from the branch cut [—1, 1] on the real 
axis and the two poles ±1/^/7. This time K(itjt) is negative 
for 2 close to 1 when S52 > 0. Moreover, the two poles tend 
to infinity when e ^ 1 so to avoid integrating too close to 
the poles, we integrate over the line 

(A9) 



2 = 1 + e^^'^y, < y < +00. 



In that case, the integrals are not simply bounded anymore 
for nearly parabolic trajectories. One can show, by making 
the previous change of variables, that 
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This expression is ill behaved for nearly parabolic trajec- 
tories, as the exponential factor diverges whereas the T- 
integral vanishes. In order to avoid numerical problems, in 
the case of nearly parabolic trajectories, we make the change 
of variables 

u = {tan ^ (All) 
The expression for the X- integral is then, for e — 1 <C 1: 
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Note that in terms of the true anomaly / we have u = 
— tan //2 and the exp r ession for the time can be found in 
iGevhng fc WestermanI (Il97ll ). Eq. (2.3.9). 
Here again we integrate along u = e"^''®j/, < y < oo, which 
cancels the ©(ii^) real part of the time and maximizes its 
positive imaginary part at infinity. Notice that for very small 
eccentricities, this is mainly a function of {u!b/v){e — 1). 



APPENDIX B: QUANTUM TREATMENT OF 
INFRARED EMISSION 

In Section [71 we computed the net angular momentum loss 
due to infrared emission using classical electrodynamics. 
Here we reconsider the effect with a quantum calculation. 
We assume a spherically symmetric grain for simplicity, and 
neglect vibration-rotation interaction. We will recover the 
classical result in the limit J ^ 1, which is applicable to the 
dust grains considered in this paper. 

The Hilbert space of the grain is characterized by the 
vibrational quantum numbers (generically denoted v) and 
the three rotation quantum numbers J, K, M, where K is 
the projection of angular momentum onto the grain ^-axis. 
The energy levels are given by 



E 



J,K,M,Ti 



21 



(Bl) 



where / is the grain moment of inertia and E'^ is the vibra- 
tional energy. The rotational wave functions are 



J,K,M 



(x) = 



2J + 1 , . 

g^2 ^K,M\X)i 



(B2) 



where x ~ {^:4':^) G SO(3) is the set of Euler angles, 
Stt^ is the volume of SO(3), and D'' is the passive rota- 
tion matrix in the spin- J representation, i.e. D'lf^ j^j^ = 

(</ A'll grain |JM2 1ab). 

Spontaneous infrared vibrational transitions are possi- 
ble from vibrational state v to v'; their rate is given by 

AjK,M.v-^J'K',M'-v' ~ 



4:{Ej,K,K 



J',K',M'.v' 



x\{j',K\ M',v'\fj.\J, K, M,v)\'' ,(B3) 

where fi is the electric dipole moment operator. In the ab- 
sence of vibration-rotation interaction, we may take the op- 
erator fi to depend only on the vibrational quantum num- 
bers and on the rotation matix R(x) that converts grain- 
fixed to lab-fixed coordinates: 

{j\K\M',v'\n\J,K, M,v) 

= {j',K\M'\R{x)\J,K,M){v'\tJ.<^'^\v). (B4) 

Here /x'^' is the dipole moment in grain coordinates. 

The transition rates can be determined by writing /x'^^ 
in the polar basis. 



and fij^l 



V2 



in which {/im }m=-i transform in the L = 1 representation 
of SO(3). Written in this basis, the rotation matrix R(x) is 
the inverse of D^(x), which for unitary is the same as 
the Hermitian conjugate: 

{J',K',M',v'\^im\J,K,M,v) 



{.;', K', M'\Dl:,^^{x)\J, K, M){v'\fil^}\v). 



(B6) 



The first matrix element can be evaluated by the three ro- 
tation matrix integral, 

{j',K',M'\Dl:,^^{x)\J,K,M} 

_ V(2J'-|-1)(2J + 1) 

^ 87r2 

= V(2,/' + 1)(2J + i)(_i)^f'+™'+Af'+™ 



J' 


1 


J 


-K' 


-m' 


K 


J' 


1 


J 


-M' 


—m 


M 



(B7) 

This transforms the spontaneous decay rate, Eq. HB3|I . into 



X E E (-ir(«'iML^'i«) 

m — — 1 m' — — l 

J' I J 

-K' -m' K 



X 



J' 

-M' 



J 

M 



(B8) 



We would now like to find the net decay rates to states of 
different J'. To do this, we assume the grain is randomly 
oriented, i.e. we average over initial projections K, and sum 
over final projections K' . Using the 3j symbol orthonormal- 
ity relations, one obtains 

-(2/ -hi) 



A 



^Ej,^ - Ej, 



J,M.v~yJ'.M'.v' — 



J'.Tj' 



E|H 



J' 

-M' 



E 



1 



J 
M 



(B9) 



(The terms mixing different values of m' are eliminated by 
orthogonality relations.) The summation over m of course 
has at most one term, with m = AM = M' — M . 

We are interested in the net angular momentum loss, 
which is most easily obtained by taking an initial state with 
M = J. There are then six possible values of AJ and AM, 
constrained by selection rules (A J, AA4" — —1,0,-1-1) and 
the restriction AM ^ AJ. The branching ratios are con- 
strained by (i) the energy difference factors in Eq. (|B9p . (ii) 
the factor of 2J' -I- 1, and (iii) the 3j symbol. We consider 
each. 

The energy factors do not depend on AM. If we take 
natural frequency u = {E'i — E'^,)/h, then the energies dif- 
ferences are given by 



Ej.v — Eji^t = hv - 



|j(2J-Kl-H AJ)AJ. 



(BIO) 



The classical grain rotation rate is a; = hJ/I. In the limit 
of J 3> 1 and Lu <ti v, the energy difference is proportional 
to 1 — a;AJ/27r!/, so the cube of the energy difference is 
proportional to 1 — Su!AJ/2nv. 
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The square of the 3j symbol, multiphed by 2 J' + 1, can 
be directly evaluated for the six cases of interest. It is: 

2J~1 
2J + 1 
1 



J+1 

J 
J+1 

1 



(J+1)(2J + 1) 
1 





-1, AM 


= -1, 


AJ = 


0, AM = 




AJ = 


0, AA/ = 


0, 


AJ = 


+1, AM 


= -1, 


AJ = 


+1, AM 


— 0, and 


AJ = 


+1, AM 


= +1. 



J + 1 

1 AJ = +1, AA/ = +1. (Bll) 

By multiplying these relative probabilities by 1 — 3ti;A J/27r!/, 
it is easily seen that the average (AM) is exactly zero if 
Lo = 0. Therefore the leading contribution to (AM) can 
be obtained by taking the large-J limit of the 3j symbols. 
Transitions with AJ ^ AM are suppressed by powers of J 
in Eq. (IB11|I . so one has three available transitions: AA^ = 
-1,0, +1, AJ = AAf. Since the factors in Eq. (|BTT|) go 
to unity, the branching ratio for these three transitions is 
determined entirely by the energy factor: 

P(AM) = i - 7:^AM. (B12) 

3 ZTTU 

This implies an average loss of 2-component of angular mo- 
mentum 

(AA/)=- — . (B13) 

In particular, we may find the ratio of angular momentum 
loss to energy loss (hu), which is 

-J -^r 2^- 

With the normalization of Eq. (I150|l and this ratio, one re- 
covers Eq. pSljl . 
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